Quantum tele-amplification with a continuous variable superposition state 
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Quantum superpositions of coherent states, often referred to as Schrodinger cat states, play an 
important role in fundamental tests of quantum theory and in quantum information processing 
and communication tasks, such as coherent-state quantum computing (CSQC) [IHSL metrology 
[6] [7] , and quantum repeater [5]. Optical cat states are now routinely generated in the laboratories. 
However, real use cases of these tasks have not yet been demonstrated. Here we demonstrate a basic 
CSQC protocol, where a cat state is used as an entanglement resource for teleporting a coherent 
state with an amplitude gain, extendable to coherent state qubits. We also show how this can be 
extended to a loss-tolerant quantum relay of multi-ary phase-shift keyed coherent states. These 
protocols could be useful both in optical and quantum communications. 



Among the various optical implementations of quan- 
tum information processing (QIP), coherent-state quan- 
tum computing (CSQC) is of special interest for enhanc- 
ing the performance of advanced optical communications, 
where information is encoded into coherent states. These 
are the only pure states that can be transmitted intact 
through a lossy channel since they are eigenstates of the 
annihilation operator, a\a) = a\a). Hence, a simple clas- 
sical encoding with coherent states - that is, without 
quantum entanglement - can be the optimal strategy of 
the transmitter to achieve the ultimate capacity of a lossy 
optical channel [9j. On the receiver side, the sequence of 
coherent-state pulses should be decoded by projecting 
them onto a set of entangled coherent-state bases so as 
to induce quantum interference in the measurement pro- 
cesses. This can be implemented by first transforming, by 
CSQC, the received codeword into appropriate superpo- 
sitions of several possible codewords, followed by a con- 
ventional separable measurement [10 j. This scheme can 
realize communication with larger capacity for a given 
input power, beating the conventional limit of optical 
communications with homodyne receivers [11] . Although 
practical implementation of CSQC remains a big chal- 
lenge, recent progress in generation [T2HT6] and manipu- 
lation [TTl [TSj of optical Schrodinger cat states makes it 
realistic to implement its basic building blocks. In this 
paper, we propose and demonstrate the first operational 
application of cat states for QIP, where a cat state is used 
as the entanglement resource for teleporting a coherent 
state with an amplitude gain. We also propose a new 
application of this to quantum key distribution (QKD), 
namely a loss-tolerant quantum relay of multi-ary phase- 
shift keyed (M-PSK) coherent states that does not as- 
sume a trusted node. We present its proof-of-principle 
demonstration with binary PSK states. 

The basic scheme of teleportation from Alice to Bob of 



a cat state qubit \i/j) a = c+ \a) A + c_ \—a) A , which is a 
variation of the schemes in [HOES], is depicted in Fig. HJi. 
Bob prepares a cat state \$_(/3)) B = Af- b ~ \ ~P)b) 
(with normalization M-) and splits it into an entangled 
cat state over paths B and C via a beam-splitter (BS) 
Vbc with reflectivity Rb- He sends one part of it to Alice 
at port C. She then combines it on an R A reflectivity BS 
with her input \ip) A at port A as 

Wabc = Vac\iP) a Vbc\*) b \0) c - (1) 

She finally measures modes A and C by single photon 
detectors. By conditioning port B on her measurement 
result, Bob can restore Alice's input. 

In order to do so, we impose a condi tion on the ampli- 
tude of the resource cat state, f3 = otyj (1 — R A ) / 'R a Rb, 
such that the components at port A turn into either the 
vacuum or a no n- vacuum state. Then, when Alice's de- 
tectors register a single photon at port A and nothing at 
port C - denoted (1,0) - Bob unambiguously obtains the 
state (see Appendix |A|) 

ac(1,0|*)abc oc c + \-ga) B + c_ \ga) B , (2) 

where g = y(l — R A )(1 — Rb) / RaRb is the gain pa- 
rameter. By a simple 7r-phase shift, it can be trans- 
formed to Alice's input \ip) A: but with a modified am- 
plitude cJ = ga. This process, previously suggested in 
[3 , we will call tele-amplification. 

Unfortunately this tele-amplification is vulnerable to 
losses. Suppose the channel between Alice and Bob is 
subject to a linear loss Re- The amplitude of the resource 
cat state should then be chosen as 

P = ^R A W-R E ) a - (3) 
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FIG. 1: Scheme of quantum tele-amplification and 
quantum relay, a Tele-amplification of binary cat-state in 
an ideal lossless channel. Ra and Rb are the reflectivities of 
the BSs. b Loss tolerant quantum relay Re is the reflectivity 
of the BS which models the lossy channel, c Alice's four-port 
measurement for the case of 4PSK states. 



After conditioning on Alice's detection event (1,0), Bob's 
state gets entangled with an external mode E as 

\°)e ^ c + \-9<*) B \e) E + c - \9 a ) B \~ £ )e ( 4 ) 
with the modified gain including the loss rate Re 



I (1-R a )(1-R b ) 
R a R b (1-R e ) 



(5) 



Here e = a/(1 - R a )R e /Ra(1 ~ Re)ol. Thus, the out- 
put at Bob is generally in a decohered state. 

One can, however, see that if Alice's inputs are re- 
stricted to classical components, namely \a) or |— a), as 
in Fig. [TJd, the output state can be completely disentan- 



gled from the external mode. This means that the coher- 
ent states can be tele- amplified faithfully to the target 
states even through the lossy channel as 



\±a) A ^ \±ga) 



B • 



(6) 



This is referred to as loss-tolerant quantum relay. In this 
context, Bob plays the role of an intermediate node, re- 
stores the target states \±ga) B , and sends them to the 
terminal node, Charlie. 

This simplest binary case can be extended into M- 
PSK coherent states. Let us show it for the 4-PSK case, 
Wm) 1 ( a m = i m a, m = 0, 1, 2, 3). Bob should prepare a 
4-component cat state 



(7) 



k=0 



as a resource state, which is quite challenging, though. 
This state is beam- split, and is shared with Alice. We 
set Ra = 0.5. As in Fig. []J;, Alice performs a four-port 
single photon detection at paths A, A', C, and C on this 
state. Depending on the set of results at the four ports, 
(A, A', C, C), the inputs are tele- amplified as 



Wm) H> \ga ml 1 

m I 1 

\a m ) >-> \-ga m ) , 
\a m ) !-> \-iga m ) 



for (0,1,1,1), 
for (1,0,1,1), 
for (1,1,0,1), 
for (1,1,1,0). 



(8) 



Thus, the simple tele-amplification is performed for the 
result (0,1,1,1). Moreover the output state can be 
switched to another element by choosing an appropriate 
click pattern at Alice (see Appendix |B| . 

The relay function itself can also be realized in a clas- 
sical way, where Bob at the intermediate node performs 
an unambiguous state discrimination (USD) on the sig- 
nal state, reproduces an amplified state for his confident 
result, and finally resends it to Charlie. This classical 
relay cannot, however, be applied to a QKD relay node 
without the trusted node assumption. Consider QKD 
protocols which use M-ary coherent states in their imple- 
mentations, such as B92 and BB84, possibly with strong 
reference pulses [24-26 . When the intermediate node 
(Bob) cannot be trusted, being at hand of an eavesdrop- 
per (Eve), a significant portion of the information shared 
by Alice and Charlie leaks to Eve by the USD, degrad- 
ing the QKD distance and the key rate. In contrast, 
a quantum relay can be carried out in the fully quan- 
tum domain, without Bob's knowing the signal state it- 
self, though at the expense of preparing the entangled 
cat-state. Similar ideas with entangled photon pairs for 
single-photon QKD were presented in [2TJ [22] . So if the 
QKD protocol is accompanied by an appropriate entan- 
glement verification session, the scheme can be used for 
a relay without the trusted node assumption. However, 
for precise quantification of relayed QKD performance of 
the current scheme, further study on validation of the 
security proof with added reference pulses as in 
is needed. 




FIG. 2: Measurement results for the twelve tele- amplifications. The blue straight lines represent a' — got. The open 
circles represent the sets of a and gt g in Table I. The fidelities between the measured states and the targeted \gtg&) are indicated 
by the numbers in brackets. Filled circles represent the amplitudes a (a) of the coherent states that have the maximum fidelity 
with the measured input (output) states. Insets: Output state Wigner functions are shown in blue contour plots. The red 
solid and dashed circles are for the input and ideal targeted states, respectively. 



We carried out an experimental demonstration of the 
tele-amplification in the simplest case of binary PSK as 
in Fig.[T}3 to realize Eq. The resource cat state |<£_) 
was generated by photon-subtraction from a squeezed 
vacuum with anti-squeezing along the real axis in phase 
space (see Methods). Bob's BS was set to Rb = 0.1. For 
a given desired gain g, we varied Ra according to Eq. 

The resource cat-state amplitude /?, experimentally 
tuned by the squeezing level, was then set by Eq. 
The detector at port C was omitted with negligible effect 
on the outcomes since the events (A,C) = (1,1) would be 
rare for our experimental setup. Bob's output state was 
characterized by homodyne tomography. 

We tested in total twelve settings as summarized in 
Table I. The five different input amplitudes a were real 
and ranged between 0.35 and 1.4. The protocol was 
carried out only for \a) because the outcome for |— a) 
would be trivially identical. The measured results are 
shown in Fig. |2j The blue straight lines are gain curves 
a' = ga in the (a, a') diagram. The open circles plot- 
ted along these lines represent sets of a and g tg in Table 
I. The Wigner functions of the tomographically recon- 
structed tele-amplified output states p out are shown as 
blue contour plots in the insets. One contour level is 
highlighted for comparison with the targeted states |#tg<^) 



(red dashed) and the actual input states p- in ~ |^)(^| ( re d 
solid, also characterized by homodyne tomography). The 
discrepancies between p out and |#tg<^) are due to imper- 
fections, including the deviation of the photon-subtracted 
state from the ideal resource cat, losses, impurity, and 
Alice's use of an on/off detector instead of two single- 
photon detectors. For each setting, we calculate which 
coherent states \a) , \a') have the highest fidelity with 
the measured input and output states, respectively, that 
is, a = argmax 7 (7|/5 in |7) and a! = argmax y (Y \p ou tW). 
These (a, a') pairs are marked as filled circles. 

Despite the imperfections, the tele-amplification suc- 
ceeded with high fidelities (gtg<^|Pout|#tg<^) between 0.89 
and 0.95 as shown next to each inset. The obtained 
amplitudes (filled circles) are close to the targeted ones 
(open circles) in almost all cases. The success proba- 
bilities were in the range 0.3% to 0.65% (see Methods). 
For larger a', the Wigner function shapes are slightly 
elongated due to the larger squeezing needed to produce 
those states. We note that our experimental settings were 
not fully optimized by taking into account spectral mode 
mismatch between the resource cat and input coherent 
states as well as between the APD and homodyne detec- 
tors. Had we done it, we estimate the achieved fidelities 
to have been 0.94-0.99. 
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Ra 


Re 


1 


0.35 


3.0 


1.11 


0.50 





2 


0.35 


2.2 


0.81 


0.65 





3 


0.50 


2.2 


1.16 


0.65 





4 


0.50 


1.5 


0.79 


0.80 





5 


0.71 


1.5 


1.12 


0.80 





6 


0.71 


1.0 


0.75 


0.90 





7 


1.00 


1.0 


1.05 


0.90 





8 


1.00 


0.76 


0.80 


0.94 





9 


1.41 


0.76 


1.13 


0.94 





10 


1.41 


0.50 


0.74 


0.97 





11 


0.35 


3.0 


1.11 


0.50 


0.8 


12 


0.35 


3.0 


1.11 


0.83 


0.8 



TABLE I: Desired tele-amplification for the twelve settings of 
input coherent states and gains. g tg is the targeted gain. 



The settings #11 and 12 had an additional 80% loss 
(Re = 0.8) in the channel from Bob to Alice. In #11, 
Ra = 0.5 for the original lossless setting (as in #1), while 
in #12, Ra = 0.83 as optimized according to Eq. ([5]), re- 
sulting in success probabilities of 0.17% and 0.11%T The 
fidelities with the target state are as high as 0.839 and 
0.872, respectively, as compared with 0.901 in the loss- 
less case. This demonstrates the loss tolerance of the 
protocol. 

Teleportation of a cat state qubit as in Eqs. ([!} 
A5) is a prerequisite for CSQC. Interestingly, the tele- 
amplification allows to convert between different ampli- 
tude qubit bases. Although we previously generated such 
arbitrary cat qubits [T7J [23] , it was not feasible to tele- 
amplify them with the current setup since three simulta- 
neous APD clicks would be needed. Instead we simulated 
this protocol by accurately modelling the current exper- 
iment including all relevant practical imperfections (for 
details, see Appendix |f|) . A cat-state qubit in the Bloch 
sphere is 



\^(a,0,(/>)) = c + \a) 



-a) 



= cos-|$+(a)> + e^sin-|$_(a)>, 

where \<f>±(a)) = Af±(\a) ± |— a)) are the even/odd cat 
states with M± = l/\/2(l ± e~ 2a2 ) and c± = A/"+cosf ± 

For a given input state, \ip(a, 0, 0)), the model returns 
an output state p a ,e,(i>, and we calculate the average fi- 
delity of the teleported state with the target state: 



TT av g _ / 

a ^ a ' J 



= / d(t)d0 



sin# 

47T 



(/))\p a ,e,cf> |^(a',0,0)) . 



In Fig. [3] we plot these average fidelities as a function of 
input amplitude a and output amplitudes a' . We see that 
for a wide range of parameters it is possible to surpass 
the classical teleportation limit of 2/3. 




FIG. 3: Simulated average qubit teleportation fidelity. 

Average fidelities as a function of input (a) and output state 
(a) amplitudes for simulated qubit state tele-amplification. 
The red curve labelled "F = 2/3" indicates the classical tele- 
portation bound. 



Finally we compare our scheme (Fig. [T^i) with the 
quantum scissor unit used in quantum noiseless ampli- 
fiers with single-photon ancilla [27-29 j. The schematic 
difference is only whether the ancilla is a single-photon 
state or a cat state. The main functional difference is the 
fidelity to the target state. In the noiseless amplifiers, 
the fidelity and the success probability rapidly decrease 
as the amplitude a increases, because the output state 
stays in the truncated space spanned by the vacuum and 
single-photon states. Hence the allowed region for a is 
limited for reasonable success rates. Our scheme attains 
in principle perfect fidelity for arbitrary amplitudes with 
higher success rates, though at the expense of restricting 
the input states to finite elements of coherent states. 

In summary, we presented tele- amplification of coher- 
ent states and a loss-tolerant quantum relay as the first 
operational application of optical cat states. The scheme 
is an essential building block for CSQC as well as quan- 
tum communications. 



Methods 

We generated the squeezed vacua at 860 nm wave- 
length from an OPO (optical parametric oscillator) con- 
tinuously pumped with pump parameters between 0.15 
and 0.31, corresponding to f3 values of 0.78 to 1.15. We 
tapped off 5% of the squeezed beam on a BS and guided 
it to an APD. A click of the APD heralded the subtrac- 
tion of a photon from the main beam [I4j [23] . The state 
thus generated is a close approximation to the odd cat 
state 1$-), and has been shown to provide near-perfect 
teleportation performance [30 . 

Whenever Alice's APD clicked simultaneously with the 
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heralding signal of the single-photon subtraction for the 
resource cat-state generation, the tele-amplification was 
successful, and we recorded a trace of the homodyne sig- 
nal of Bob's output state. The success probability is 
given by the ratio of the simultaneous click rate (~3-28 
to the photon subtraction click rate (^1000-4500 
It is mainly limited by detector and spectral fil- 
tering efficiency. To build the homodyne tomogram, we 
repeated this procedure 6000-24000 times for each fixed 
input state, with the local oscillator of the homodyne de- 
tector locked at phases -180°, -150°, . . . , 150° with re- 
spect to the input state. Note that the protocol succeeds 
as a single shot for an unknown input state - the repeated 
measurements with identical inputs are only needed for 
characterizing the process by homodyne tomography. 

Alice's input states were independently characterized 
by homodyne tomography at port C by setting Ra = 1. 
To determine the input states accurately just at Alice's 
BS, we correct their reconstruction for the detection ef- 
ficiency and the propagation losses from that point to 



the homodyne detector. This total efficiency amounts 
to 88%. Likewise, in the reconstruction of Bob's output 
states we correct for the overall detection efficiency of 
94% but not for any propagation losses. 

A more detailed description of the experimental setup 
and the state characterization can be found in Appendix 

m 
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Appendix A: Tele-amplification of a binary through a lossless channel. The whole three-mode state 

component cat-state after the beam-splitting operation in Fig. [T^i is given by 

We first describe the tele-amplification of a binary 
component cat state \ip(a)) A = cq \ol) a + c\ |— a) A 



l*W = VacM<*)) a Vbc\*)b\0)c 
= Afc 



\fl ~ R A a - ^RaRb/3) \A - RbP 



R A a - - R a )Rb/3 



c 



Afc a/1 " R A a + v'RaRbP) -y/\-R B P 



R A a + V(l - R a )RbP 



c 



- a/1 - R A a - ^RaRbP) a/1 - RbP 

I A IB 



R A a - v/(l - R a )RbP 



c 



x/1 - Rao + VR A Rb/3) -y/l-R B P) + - R A )R B P 



c 



r 



(Al) 



where Af = 1/ \/2[l — exp(— 2|/3| 2 )] is the normalization 
of the resource cat state \<&) B . We now impose a condi- 
tion on the amplitude of the resource cat-state 



/3 



l-RA t 
RaRb 



(A2) 



such that the components at port A turn into either of 
the vacuum or non-vacuum states as 



I*) 



ABC 



Mc \0) A \ga) B 
AT c 2^1 -R A a 



1 



-2y/l - R A a 



a 

'Ra i c 

A \-9 a ) B 
\9<*)b 



1-2R A 
— a 
VRa I c 

1-2R A 



ATct \0) A \-ga) 2 



a 

'Ra I c 



with the gain 



g = ^ (\ - R A )(± - RbMRaRb 



'Ra I c 
(A3) 

(A4) 



Alice then performs single photon detection on paths 
A and C as shown in Fig. [T^i, and selects single photon at 
port A and nothing at port C - denoted (1,0). Then Bob 
can unam biguously exclude the first and fourth terms in 
Eq. ( A3 ) , and has the state 



ac(1,0|*)abc \4>(-ga)) 



B ' 



(A5) 



In the case where the channel between Alice and Bob is 
subject to a linear loss with the rate Re, one can consider 
an external mode E. Bob chooses the cat-state amplitude 
as 



1-R A 



R a Rb(1-Re) 



(A6) 



The whole state before Alice's measurement is 
Wabce 

= VacW a VecVbc\$) b \0)c\0)e 
1 



= Afc \0) A \ga) 1 



^r A a /c h£)E 



■Nco 2y/l - R A a) \-ga) 



l-2R A 



+J\fd -2y/l - R A a) \ga) B 

I A 



-Afci|0) A \-ga) B 
with the gain 



VRa i c 
l-2R A 



fRl 



-a 



l (1 - R A )(1 - R B ) 



(A7) 



(A8) 



y R a Rb(1 — Re) 
and e = - R A )R E /R A (l ~ Re)ol. 

Appendix B: Extension to multi-ary coherent states 

The binary case can be extended to M-ary phase-shift- 
keyed coherent states \a m ), where 



a rn = au , u = e 



2iri/M 



(Bl) 



Here ao = a is taken to be real. The states are generated 

by 

\a m ) = V m \a ), l/ = exp(^n). (B2) 
An input state at Alice is 

M-l 

\^)a = Yl ° m \ a rn) A - (B3) 
ra=0 
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Bob prepares a cat state for the entanglement resource, 

M-l 



m B = E b k\f3 k ) B - 



(B4) 



k=0 



In order to see what kind of cat state is appropriate, we 
introduce the orthonormal basis {|cj m )}, which diagonal- 
ize the generating operator as 



M-l 

V=J2 u k \oo k )(uk\ 

k=0 



Then 



M-l 



(B5) 



(B6) 



m=0 

Using the orthogonality property 

M-l 



^ m ~ n)k = MS min , 



(B7) 



k=0 



we have 



M-l 

The normalization condition {uj m \uJm) = 1 means 

M-l M-l 

I MP) I 2 = W E E u(fc '" fe)m • (B9) 

/c'=0 A:=0 

Because the matrix {{Pw \Pk)) is cyclic, 

M-l 

M 



M-l 

I" = E U ~ km (A)|/3fc>. 



(BIO) 



/c=0 

Define a quantity 

A m = M|( Wm |/?)| : 

Then one can see 

M-l M-l 



^77 



(Bll) 



(B12) 



m=0 



ra=0 



These eigenvalues are related to the mean photon number 
of the basis states as 



(u) m \n\uj m ) - — . 



(B13) 



To maximize the success probability of Alice's measure- 
ment, one should use the \uo m ) which has the maximum 
photon number, namely the minimum eigenvalue, for the 
entanglement resource. For relatively smaller |/3|, it is 
the \uj M -i)- 



Now let us see the case of M = 4 (u = i). The cat 
state for the entanglement resource is 



M-l 



l*>B = \<**) B = ^E^ \h 



The eigenvalues are given by 



(B14) 



k=0 



A = 2e~P 2 (cosh j3 2 + cos j3 2 ) 
Ai = 2e-^ 2 (sinh/3 2 + sin/3 2 ), 
A 2 = 2e~^ (cosh /3 2 -cos /3 2 ) 



(B15) 



A, 



2e-^ (sinh/3 2 - sin/3 2 ). 



The above state is beam- split into paths B and C, the 
component of mode C is sent to Alice through a lossy 
channel, and then combined with the input at path A. 



Bob chooses the cat-state amplitude as Eq. ( A6 ) . The 



whole state before the measurement is given by 

3 3 



1*1 



BACE 



ra=0 £:=0 



\y/l-R A a{u m -u k ) 



a 



\eu k ) 



Ra 

E 



R A u m + (1 - R A )u k 



c 

(B16) 



with the gain given by Eq. ( A8 ) 



We set Ra = 0.5. Alice further introduces additional 
modes A' and C to implement the four-port single pho- 
ton detection as shown in Fig. [TJd. The state before the 
detection is given by 



I*) 



BAA'CC'E 



m=0 k=0 
u rn — U k 



,m _ „.ra+l 



+ u k + u k+1 



2V2 

m+1 



+ u k - u k+1 



2V2 



VI -Re 



c> 

(B17) 



When th e loss can be neglected (Re = 0), the input 
of Eq. (B18) can be faithfully tele- amplified to the 

target state 



M-l 



\9^)b = Yj Crn \d a m) 



B ' 



(Bis) 



rn=0 
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by selecting a set of Alice's measurement result as (A, 
A', C, C')=(0, 1,1,1), namely no count at port A while 
single photon counts at port A', C, and C\ 

In the lossy case, it is impossible to teleport a super- 
position state faithfully. However, when an input is re- 
stricted to a classical state drawn from the set \a m ), then 
the tele-amplification to the target pure state is possible. 
Actually depending on a set of the results at the four 
ports, (A, A', C, C), the inputs are tele-amplified as 



\ot m ) ^ \gOLm) , for (0, 1, 1, 1), 

Wm) *-> \igot m ) , for (1,0, 1, 1), 

Wm) i-> \-goL m ), for (1,1,0, 1), 

\a m ) \-> \-igam) , for (1,1,1,0). 



Appendix C: On/ofF detection at Alice 



(B19) 



II off = |0> <0| and n on = I — |0) (0|. Then the tele- 
amplification described in the previo us se ction should be 
corrected slightly. For example, Eq. (A5) for the binary 
case becomes 

ABC (m±fl^\y)ABC = m-ga)) B (^(-ga)\ 
+ tanh(|)|^(-^a)> B (^(-^a)| (CI) 



where \ijj(—ga)) = cq \a) — c\ \ and a = 2^1 — Ra®- 
The second term is the correction. When a is small, the 
coefficient of the second term is small as tanh(a/2) ~ 
a/2. In this regime, the tele- amplification would approx- 
imately work with on/off detection. However, in general, 
the second term cannot be ignored. 

If the input state was a coherent state, i.e., cq = 
or ci =0, the above state would become a pure coher- 
ent state with the gain g, and on/off detection would be 
sufficient. 



In our experiment, Alice's measurement is imple- 
mented by avalanche photo diodes (APDs) instead of 
ideal "single photon detectors" that discriminate be- 
tween "0", "1" and "2 or more" photons. APDs can- 
not, however, dicsriminate photon numbers, but dis- 
tinguish merely the vacuum or non-vacuum state, i.e. 
"off' or "on". They are represented by the operators 



Appendix D: Success probability 

1. Loss tolerant quantum relay 

In the binary case, the input is either of \a) and \—a). 
The whole state before Alice's measurement is 



|*±) = V AC \±a) A V EC V BC \<$>) B \0) c \0) E 

1 



= ±M\Q) A \±ga) 1 



±2y/l -R A a) |TWi 



A-2R A 

± — a 



iRa 



\±a') E (Dl) 



r 



The success probability of the tele-amplification In the case of 4-PSK states, the state before Alice's 
|±a) A i-> \±ga) B is given by the expectation value of measurement is given by 

ftio = n4 ® n& as 



4 2 l- amp = ^(* + |niol* + ) + ^(*-|niol*- 



exp 


(1-2R A ) 2 2 
Ra a _ 


— exp 


Ra 


2 (l - exp 


2(1-Ha) a 2 
RaRb(I-Re) 


) 



.(D2) 



k=0 



U m +U K 



u m _ u m+l + u k + u k+l 



2V2 



u m + u m+l + u k _ u k+l 



2V2 



Re k 



VI -Re 



(D3) 



The success probability of \oi m ) A i-> \ga m ) B is given by 
the expectation value of 



ftom = iioff ® n£ ® tf on ® tiZ 



(D4) 
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FIG. 4: Success probabilities for the case of the BPSK co- 
herent states. The solid lines are for the quantum relay with 
several cases of the amplification gains, while the dashed line 
is for the measure-resend strategy with the unambiguous state 
discrimination. 
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FIG. 5: Success probabilities for the case of the 4PSK co- 
herent states. 



as 



d(4) 

Tele — amp 



^<* ro |noni|* ro ) 

k=0 

(*o|n in |* ) 

(l_ e -« 2 /2)2( 1 _ e -a 2 ) 



4A 3 ( 



Rb (1 — Re) 



where 



Xs(x) = 2e x (sinh x — sin x) . 



(D5) 



(D6) 



He tries to discriminate them unambiguosly without er- 
rors, but at a finite success rate, and then prepare a 
target amplified state \gy/l — Re^tti) for the measure- 
ment result m. The success rate is well known for this 
kind of equally probable symmetric states [1 j. Denoting 
1 7m) = | Vl — REOim) and using the eigenvalues and the 
diagonalizing vectors of the density operator 



M-l 

P= Yl l 7m ) ^ 



M-l 



A « 

m=0 m=0 

the success rate is given by 



min Xh. 

k 



PlJSD 

The signal states are represented as 



(D7) 



(D8) 



|7 m ) = \A - Reoit, 



M-l 



\ k u mk \ui k ) . 



(D9) 



fc=0 



The detection operators are given by 
A 



n »> = ^PusD\i^){ii\ 



(D10) 

for the signal state |7 m ), using the reciprocal states 



1 M-} mk 

^ kT ^ 

where A = ^2 k A^ 1 . They satisfy the orthogonality rela- 
tion 



<7ml7m') 



(D12) 



The operator for the inconclusive result is given by 

M-l 

fl F =I - Yl ft m. (D13) 

rn=0 



3. Numerical results for BPSK coherent states 

The success probabilities are shown in Fig. [4] for the 
BPSK coherent states, and in Fig. [5] for the 4PSK co- 
herent states, in comparison with that by the measure- 
resend strategy. The quantum relay can attain higher 
success rates than that by the above measure-resend 
strategy. In the case of 4PSK states, the difference is 
as high as ten times. 



2. Measure-resend strategy 

The task to relay attenuated coherent states to the re- 
ceiver, converting them faithfully to the target amplified 
states, can also be realized by a classical strategy. A 
typical one is a measure-resend strategy. In the interme- 
diate node, Bob has attenuated states — RE&m)}- 



Appendix E: Experimental details 
1. Setup 

The most relevant elements of the experimental setup 
are sketched in Fig. [6] and described in the following. 
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FIG. 6: Sketch of the experimental setup for tele- amplification of a coherent state from Alice to Bob using shared entanglement 
in the form of a photon-subtracted squeezed vacuum state. 



Resource state generation 

The output of an 860 nm continuous- wave Ti:Sapph 
laser (a) is phase- modulated by an EOM (b) for Pound- 
Drever locking of the SHG and OPO cavities. Parts of 
the beam are tapped off for use as local oscillators, coher- 
ent input beam, probe beam and cavity locking beams, 
but the main part is frequency-doubled in the second 
harmonic generator (c). The 430 nm output of this SHG 
pumps a bow-tie configuration optical parametric oscil- 
lator (OPO, d) with a PPKTP crystal and a HWHM 
bandwidth of 7/27T = 4.5 MHz. The down-converted 
light leaving the cavity is in a squeezed vacuum state 
(e). A probe beam (f) is injected into the OPO for the 
purpose of locking phases and filtering cavities further 
downstream. 

The whole experiment is running in an alternating 
lock/measure cycle at a rate of 10 kHz. The probe beam 
is thus switched on and off by double-passes through two 
acousto-optical modulators (only one shown in the figure) 
that are driven for only 20% of the 10 kHz cycle time, as 
indicated by the blue pulse diagram (g) (this diagram is 
repeated for other relevant parts of the setup). 

In order to act as a phase-reference for the squeez- 
ing, the probe beam is locked in phase with the squeezed 
quadrature by observing its classical parametric am- 
plification in the OPO through a <1% tap-off of the 
OPO output (h). To obtain an error-signal, the probe's 



phase is dithered on a piezo-mounted mirror by a micro- 
controller unit (Arduino) which also processes the de- 
tected signal and provides feedback to lock the phase (i) . 

A photon is subtracted from the squeezed vacuum at 
random times by the detection on an avalanche photo- 
diode (APD) detector, placed after a 5% tapping beam- 
splitter and two frequency- filtering cavities (j). The APD 
is protected from the strong probe beam by an AOM that 
directs the OPO output to the APD only during the in- 
tervals when the probe beam is switched off. The result- 
ing photon-subtracted squeezed vacuum (PSSV) state (k) 
is the resource of entanglement in our protocol, after it 
is split into two modes propagating towards the Alice 
and Bob sections of the setup. In the description of the 
protocol in the main section, a fraction Rb — 0.1 is re- 
flected towards Alice. In the actual experiment (and in 
the setup sketch), Alice's fraction Rb is actually taken 
as the transmitted part of a variable beam-splitter fixed 
at 90% reflection (1). Bob's share of the entangled state 
is directed towards a homodyne detector for output state 
analysis. 

The amount of squeezing produced by the OPO deter- 
mines the amplitude of the cat-like PSS V state. It is reg- 
ulated by the pump parameter e = v^pump /^threshold- 
To find what cat amplitude f3 a given pump parame- 
ter corresponds to, we model the PSSV as in Ref. [23] 
of the main text. Here, we include the OPO temporal 
correlations, the temporal modes of the APD and homo- 
dyne detection (described in the following section), the 
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FIG. 7: Wigner function for one of the PSSV states used 
as approximations for the odd cat resource state, in this case 
with e = 0.20. The upper is experimentally generated and to- 
mographically reconstructed, while the lower one is obtained 
from our model. The fidelity between them is above 98%, 
showing the validity of the model. 




0.1 0.2 0.3 0.4 

OPO pump parameter e 



FIG. 8: Relation between OPO pump parameter e for the 
production of realistic PSSV states and the amplitude of the 
cat state that maximizes the mutual fidelity. The 

blue curve indicates these optimal /3 amplitudes, while the 
dashed red curve shows the corresponding fidelities. In the 
experiment, we used e in the range 0.15-0.31. 



96% escape efficiency of the OPO, the 95% propagation 
efficiency towards the beam-splitter (1), the 5% tapping 
ratio, the ~10% overall APD detection efficiency, and the 
filtering bandwidth to find a ppssv that closely emulates 
the actually produced states. Fig. [7] shows an example 
of an experimentally generated PSSV state and, for com- 
parison, the modelled state with equivalent parameters. 
We then maximize for j3 the state's fidelity with a true 
cat state, ($_(/?) | p PS sv | $-(/?)), and get the e -)> (3 cor- 
respondence plotted in Fig. [8] and thereby the f3 values 
of Table H 

Input state and teleportation 

Alice's input coherent state (m) is prepared in a con- 
figuration of two double-pass AOMs similar to that used 
for the OPO probe beam, with a strong phase locking 
beam switched on during the 20% locking part of the 10 
kHz cycle. However, instead of switching the light com- 
pletely off during the remaining 80% of the cycle, a weak 
amplitude beam is generated instead. This is done by 
switching to a lower voltage RF driving signal for the 
AOM. 

The share of the entangled PSSV state propagating 
from Bob to Alice in mode C is optionally subjected to a 
loss at a variable beam-splitter (n) before it is overlapped 
with her input state on a polarizing beam splitter (PBS) 
in orthogonal polarizations (o). A half-waveplate fol- 
lowed by another PBS (p) then interferes the two modes 
as the protocol's Ra reflectivity beam-splitter. When 
characterizing the input state, Ra is set to 1, which 
means that all of the input state is sent towards the ho- 
modyne detector (q). Otherwise, when running the tele- 
amplification protocol, Ra is set to its appropriate value, 
and the output of the beamsplitter in mode A is sent to- 
wards an APD (r) with the same frequency filtering and 
chopping configuration as that in (j). 

The detection events from the two APDs are correlated 
with digital timing electronics that pick out simultane- 
ous events and triggers the acquisition of Bob's homo- 
dyne signal at a fast digital oscilloscope (s). The detected 
photo-currents are subsequently temporally filtered on a 
PC to extract the measured quadrature values, as de- 
scribed in the following section. 



Phase locking 

For the tele-amplification protocol to work, the input 
states |=ba) should be interfered in-phase with the anti- 
squeezed quadrature of the entangled PSSV state. We 
do this by putting a normal intensity detector (t) at the 
mode C output port of the beamsplitter instead of Alice's 
homodyne detector (q). The detected interference signal 
between the probe portions of the input coherent state 
beam and the PSSV state beam is used as the input to 
an FPGA-based lock unit, which provides feedback to the 
phase of the input beam (u). When the two beams are 
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interfered at 90°, we get the desired phase relation, since 
the PSSV probe was locked to the squeezed quadrature. 

The phase of the local oscillators (LO) in the two ho- 
modyne detectors can be locked to arbitrary phases rel- 
ative to the PSSV probe by using a combination of DC 
and side-band detection of the interference between the 
LO and the probe beam. An 8 MHz phase modulation 
is applied to the LO (v), and the interference signals ob- 
served by the fast homodyne detectors (from a separate 
low-gain amplification output) are demodulated at that 
same frequency. This provides an interference signal that 
is 90° out of phase with the DC signal. In the FPGA lock 
units (u,w), the two signals are added with weighting 
factors corresponding to the desired LO phase in phase 
space, resulting in an error signal for the feedback to 
piezo-mounted mirrors in the LO beam path (in the case 
of Bob's output homodyner) or in the input beam path 
(for Alice). 

All phase locks are engaged only during the intervals 
of the 10 kHz experiment cycle in which the probe beams 
are turned on. For the remaining, the feedback signals 
are just held at their last actively set value. 



2. Quantum state tomography 

Temporal modes 

The squeezed vacuum has a bandwidth given by the 
OPO's HWHM of 7 /2tt = 4.5 MHz. Conditioned on 
an APD click at time to, the continuous- wave squeezed 
vacuum is converted into a temporally localized PSSV 
state in a temporal mode around to which, in the low- 
squeezing limit, has the form exp(— j\t — to|) El- The 
filter cavities in front of the APD, needed to remove the 
photons down- converted into the many non-degenerate 
OPO resonances, modify the temporal mode to be 



-«|t — to | 



(El) 



where nj2ix « 25 MHz is the combined bandwidth of the 
two filters, approximated by a single Lorentzian spectral 
profile. This will also be the temporal mode of the tele- 
ported output state in the low-squeezing limit. For the 
state tomography, we therefore extract a single quadra- 
ture value from the continuous photo current signal of 
the homodyne detector by integ rati ng it over a mode 
/HD,out(£) equal to the one in Eq. (El ). At higher squeez- 



ing levels, the optimal mode function is not that simple 
[3] , but in this work we stick to the simple expression for 
all squeezing levels. 

An interesting, but also complicating aspect of our cur- 
rent implementation of the tele-amplification protocol is 
that the input and output states are in rather different 
spectral modes: The input coherent state is derived di- 
rectly from the narrow-band laser, whereas the entangled 
PSSV state is in the broadband mode described above. 
At a first glance it would appear like the two modes will 



not interfere well and the teleportation will fail. How- 
ever, the spectral response of Alice's APD is very broad, 
so it is unable to distinguish the modes. Thereby it can 
be said that the detection itself induces the interference 
between the input and the entangled state. Another way 
to see it is in the time domain: compared to the cw in- 
put beam and the ~ I/7 extent of the PSSV state, the 
temporal response of the APD (~350 ps jitter) is essen- 
tially delta function-like. Within this short time window, 
almost no phase shift will occur between the different fre- 
quency components, so interference will not be destroyed. 

One problem we do get from this spectral mismatch, 
however, is the issue of which temporal mode, /HD,in(£) 
to use for the definition of the coherent input state. As 
the beam is continuous, the choice of temporal mode can 
be done arbitrarily. The photon number n- in within the 
chosen mode will be proportional to the width of the 
mode function, so to obtain a desired coherent state the 
intensity of the beam should be adjusted inversely pro- 
portional to that width. In our experiment we make the 
rather natural choice to use the /HD,out(£) mode, such 
that we observe the same temporal mode in both the 
input and the output homodyne tomography. The mea- 
sured a and a' values in Fig. [2] are therefore directly 
comparable. 

However, this mode is not the one detected by the 
APD. Since the APD with its delta function-like response 
is preceded by filtering cavities which act as delays for 
incoming fields, its temporal mode can be approximately 
described as a single-sided exponential decay, with time 
constant given by the filter bandwidth, 



/apd(*) oc e 



-K\t — to\ 



H(to - 1), 



(E2) 



with H(t) being the Heaviside step function. Because the 
PSSV entangled state is temporally localized, as opposed 
to the input state, the ratio of the photon numbers of the 
two states, ni n /npssv will be different within the differ- 
ent modes fwD,in(t) and /apd(£)- There is therefore a 
mismatch between the input state amplitude, a, that we 
expect to have and the amplitude actually seen by Alice's 
APD, which is the one to induce the teleportation. Thus, 
the /? and Ra values that we experimentally adjusted to 
match a given a were actually not optimal, and this re- 
sulted in output-to-target fidelities that were lower than 
we could have otherwise obtained. In a possible follow-up 
experiment, it would be advisable to consider this issue of 
the input state amplitude in more detail. Simulations in- 
dicate that with optimized settings, fidelities could have 
reached 0.94-0.99. 



State reconstruction 

For a given realization of the tele-amplification, we con- 
struct a homodyne tomogram of the output (and input) 
state by repeating the state preparation, on/off detec- 
tion and conditional homodyne detection multiple times, 
with the LO phase of the homodyne detector fixed at 
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FIG. 9: Example tomogram, showing the 12 x 2000 quadra- 
ture values of the homodyne measurement of the output state 
in tele-amplification #5. 



we set up a model for the protocol, using Wigner function 
formalism. 

As the input state to be teleported, we took a pure 
qubit state of the above form. The initial resource state 
was a squeezed vacuum state with appropriate squeezing 
levels. In the experiment, the squeezed vacuum state 
within the homodyne-observed mode /HD,out(£) is n °t 
pure. The impurity due to this mode selection can be 
modelled quite well by propagating the initially pure 
squeezed vacuum through a 92% transmission beam- 
splitter. The losses suffered by the photon-subtracted 
squeezed vacuum were similarly modelled by virtual 
beam-splitters, taking account of the 96% escape effi- 
ciency of the OPO, the 5% tapping ratio for the photon 
subtraction, and the 95% propagation efficiency towards 
the separating beam-splitter. Alice's detector was mod- 
elled as an on/off detector with 10% efficiency, roughly 
corresponding to our APD's detection efficiency and the 
transmission of the spectral filters. 



various angles. After filtering the oscilloscope traces 
with the chosen temporal mode, as described above, the 
obtained quadrature values are normalized by vacuum 
traces recorded under the same experimental conditions 
while we also pay attention to proper offset correction 
of the traces, which can be particularly tricky for the 
measurement of the input coherent state. That gives us 
a homodyne tomogram like the one shown in Fig. [9] 
From this we reconstruct an estimate of the underlying 
quantum state using the maximum likelihood method [4] . 
As mentioned in the Methods section, we correct for the 
non-perfect detector efficiencies in order to get the most 
accurate characterization of the protocol. 

The phase values in the figure indicate the relative 
phase between the local oscillator and the OPO-injected 
probe beam. The probe beam is locked to the squeezed 
quadrature of the PSSV state, and Alice's input coherent 
beam is locked at 90° to the probe beam. Since we define 
our phase space in such a way that the anti-squeezing is 
aligned along the x-axis and the input states have real 
amplitudes (i.e. also along the x-axis), the 90° phase of 
the LO should correspond to the ^-quadrature. We there- 
fore rotate the reconstructed quantum state by —90° in 
phase space - the free choice of global phase. 



Appendix F: Modelling of qubit teleportation 

To simulate the performance of our tele- amplifier setup 
in the case where the input is an arbitrary coherent state 
qubit 



\i>{a,e,4>)) = cos- |*+(a)> +e^sin- |$_(a)> , (Fl) 




n: 



i 



FIG. 10: Bloch sphere map of the fidelities between modelled 
and targeted outcomes of the \ip(ot,0,<fi)) —> 0,<j))) tele- 

amplification, in this case for a — 0.4 and a — 0.6. The 
averaged fidelity here is 77%. 



For a given input amplitude a and desired output am- 
plitude a' ', we simulate the tele-amplification process for 
168 evenly distributed qubit states on the (0, <p) Bloch 
sphere and calculate the fidelity between the output 
states and the targeted states \ip(a' : 0, <p)). This results 
in a "fidelity map" like the one in Fig. 10 for every (a, a') 
setting. It is clear that the teleportation works best for 
coherent state inputs (near 100% fidelity) and for states 
near the North Pole (which is the even cat) and not very 
well for states near the South Pole (odd cat). By averag- 
ing over the Bloch sphere, we obtain an average fidelity 
for the outcomes of the protocol for the given (a, a') pair, 
giving one value for the average fidelity plot in Fig. [3] 
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